Effective models of two-flavor QCD: from small towards large 
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We study effective models of chiral fields and Polyakov loop expected to describe the dynamics 
responsible for the phase structure of two-flavor QCD. We consider chiral sector described either 
using linear sigma model or Nambu-Jona-Lasinio model and study how these models, on the mean- 
field level when coupled with the Polyakov loop, behave as a function of increasing bare quark (or 
pion) mass. We find qualitatively similar behaviors for the cases of linear sigma model and Nambu- 
Jona-Lasinio model and, relating to existing lattice data, show that one cannot conclusively decide 
which or the two approximate symmetries drives the phase transitions near the physical point. 



I. INTRODUCTION 



It has been proposed that one can understand the interplay between chiral symmetry restoration and center sym- 
metry breaking and their implications on the phase transitions of QCD at finite temperature and density by treating 
chiral symmetry as almost exact, and center symmetry breaking as a consequence of chiral symmetry restoration and 
interactions [T]. This picture allows one to explain, in very transparent terms, why deconfinement and chiral symmetry 
restoration coincide in the theory with fundamental quarks [5] while for adjoint quarks the two phase transitions are 
widely separated [IT. For adjoint quarks at finite chemical potential this model framework also predicts interesting 
multicritical phenomena, see [4j. For QCD the standard viewpoint, supported by the physical spectrum, is that in 
the two-flavor case the chiral symmetry is almost exact whereas the center symmetry is badly broken. However, one 
might regard this picture too simple, since it does not take into account how rapidly the system departs from the 
limit of exact chiral symmetry. Consider cranking pion mass up from its physical value. There are basically two 
options: First, it could be that even for modest masses chiral symmetry remains a good approximation while the 
center symmetry becomes an approximate symmetry only for a very heavy pion. Alternatively, the second option is 
that the system departs from the chiral limit very fast, and chiral symmetry becomes very quickly a poor symmetry 
simultaneously allowing center symmetry to restore rapidly as a function of pion mass. This latter viewpoint has 
been advocated for in e.g. [5]. A new ingredient of our analysis is the treatment of chiral and Polyakov degrees of 
freedom within the same model setup simultaneously. Furthermore, in the spirit of our earlier work [S], we consider 
side-by-side two different realizations of the chiral sector. This will allow us to gain more insight into what extent the 
quantitative results of these models can be trusted. 

We consider linear sigma model and Nambu-Jona-Lasio model [7] both augmented with the Polyakov loop. There 
has been a large body of literature devoted on these models recently [5] |S1 IHl IIHI IHl HH US] ■ These models provide a 
simple and concrete realization for the intertwining of deconfinement and chiral symmetry restoration as explained in 
general terms in [1] . As discussed in detail in [6] , especially at finite densities the numerical output is quite sensitive 
on which realization of the chiral dynamics is used. Therefore the quantitative results should, when used to draw 
concusions about QCD dynamics, be interpreted with some care. 

In this paper, we continue the investigation started in j^. Here our aim is to study how these models respond when 
the amount of explicit chiral symmetry breaking is increased. For the case of Polyakov-linear sigma model (PLSM), 
we basically confirm the picture established in [5] and we find that Polyakov-Nambu-Jona-Lasinio (PNJL) model 
leads to qualitatively similar results. Our main conclusion is that the PLSM and PNJL models lead to qualitatively 
and to some extent quantitatively (within 10% or less) similar phase structure at and away from the physical value 
of TOtt. This observation strengthens the viability of these models as an effective description of the thermodynamics 
of two-flavor QCD at zero net-quark density. 

The paper is organized as follows: in section |lT] we briefly recall the basic definitions of the two models we consider 
and explain how the parameters are affected to take into account varying pion mass. Then in section [TTT] we present 
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our main results and in section HVl our conclusions and outlook. 

II. THE MODELS 

A. Models at the physical point = 140 MeV 

Let us first review the general definitions of the models which we consider. We start with the description in terms 
of parameters fixed to correspond to the physical values reproducing the observed masses for the pion and scalar 
resonance a. The PLSM model consists of the linear sigma model, a polyakov loop potential and an interaction 
between the two. The PNJL model is similar with the chiral part of the Lagrangian now corresponding to the NJL 
model. For the derivation of the grand potential we refer to the literature, e.g. and merely state the result here 

V 

~ Uchiial + Ug + flqq, (1) 

where 

C/chiral ^^{a^+TT^- V^ f - H<J (2) 

for the PLSM model and 

(^chiral - ^ (■i) 

for the PNJL model. The parameters in the above equations are fixed by the physical vacuum properties for LSM 
as H — fTTiri^, where /tt = 93 MeV and = 138 MeV. The coupling « 20 is determined by the tree level mass 
771^ = 2A^/^ + 777^, which is set to be 600 MeV. On the other hand, for NJL part we fix the bare quark mass to 
be TTiQ = 5.5 MeV and the coupling G = 10.08 GeV^'^. Furthermore, the constitutien mass M is related to the qq 
condensate as M = mo — G{qq). 

The Polyakov loop is included to both models through the mean field potential 

u, ^ u{i,t,T)/T^ = -^Kp - ^j{f+t') + (4) 

62(r) = ao + ai^+a2(^f + a3(^)^ (5) 

and the constants ai,bi are fixed to reproduce pure gauge theory thermodynamics with phase transition at Tq — 270 
MeV; We adopt the values determined in [£ , and shown for completeness in table |l] Here i is the gauge invariant 
Polyakov loop in the fundamental representation, and while one could also include other loop degrees, say, adjoint or 
the sextet, we consider here only the mean field potential of the fundamental loop parametrized to describe the pure 
gauge thermodynamics. Finally there are the interactions between the Polyakov loop and chiral degrees of freedom 
given by the potential 



where 



f]„-„ = -2NfT 



(ln(l + 3(£ + r e-(^-^)/^)e-(^-^)/^ + e-3(^^-'')/^) 



(27r)3 

+ ln(l + 3(r + ee~^E+^.)/Ty-(E+^.)/T ^ ^-^{E+^.)/T^\^ 
In the PNJL model the interaction potential ([6]) has also an additional vacuum term 

which is controlled by the cut-off A. In the above equations we have E — \/p^ + for the PNJL model and similarly 
+ g'^cr'^ for the PLSM model. In the latter case the coupling constant g is fixed to 3.3 corresponding to the 
baryon mass ^ 1 GeV. Thermodynamics is now determined by solving the equations of motion for the mean fields, 

da-'^ m-'^ df-'^ 
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^ 1.77 ■ 10-^ GeV^ 
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651 MeV 
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6.75 


-1.95 


2.625 
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-7.44 


0.75 


7.5 



TABLE I: The parameters used for the effective potential 



and then the pressure is given by evaluating the potential on the minimum: p = — fl{T,fi). The above equations are 
valid at nonzero densities, but in this work we will not consider finite chemical potential but set fj, = 0. This results 
in the additional simplification £ ^ i* . 



B. PLSM model away from the physical point 



In this work we want to explore the region in the liiodels where the pion mass differs from its observed value. For 
this we need a consistent way of setting the model parameters in the non-physical region. In the PLSM model this 
can be achieved by finding relations /Tr(?Ti^), m^{mTr) and g^mT^) consistent with the known values at the physical 
point and leaving ttItt as the only tunable parameter. In this section we present a way to do this relying on lattice 
data. 

First we want to relate the pion mass and the bare quark mass m^. From jl4j we get the relation 

mla'^ = Ai{mqa)^ + Bimga)"^, (9) 

for the pion mass given in terms of the lattice spacing a and m^. In ([9| the parameters Ai, B and 5 are fitted from 
lattice data and we use the fit made in [14] 

5^ 0.16413 

Ai = 0.82725 (10) 
1.88687. 



The lattice spacing a is determined from the equation 

/^a = 0.06672 + 0.221820 x (m,a) (11) 

also fitted in [3] to lattice data. Now to retain the relation that /jr = 93 MeV corresponds to m^r = 138 MeV we 
shift the value of /tt by a constant C = 1.18 MeV since we don't want to alter the fit parameters or the lattice spacing 
obtained in [14]. This gives us the relation 

rriqa = ((/^ -|- C)a - 0.06672) /0.221820 (12) 



with the lattice spacing a = 0.505306 GeV^ . Combining (|9| and (12 1 gives us the pion mass in relation with the 
pion decay constant. 

The sigma mass is determined through the relation 

m„^^ml + D (13) 

based on n~5| from which we also obtain the slope ^ = 0.00183 MeV^^. Requiring that TOtt = 138 MeV corresponds 
to rria = 600 MeV gives us the value of the constant D = 565.15 MeV. 

Finally, for the coupling constant g we had the relation g/^ — M^/S, where the nuclcon mass Mjv can be 
parametrized in the form 

Mjv = Mo + 4Ciml (14) 

as discussed in [16 where it has also been shown that for a good description of nucleon mass one should include terms 
up to and including O(m^) which are obtained from chiral perturbation theory at 0{p'^). We here truncate this fit to 
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0{m'^), which gives a fairly poor description of Mjv(TO7r) but the resulting effect on g is small, at 30% level and we 
will comment on the possible effect this has on our numerical results. Also, we want to keep the number of parameters 
to the minimum, so this approximation seems a re ason able starting point. To retain the relation (m^ = 138 MeV, 
= 93 MeV) <^ g = 3.3 we set the parameters in ^ to Mq = 860 MeV and Ci = 0.8 GeV"\ 



C. PNJL model away from the physical point 

In the PLSM model we first used the lattice relation ^ to connect the pion mass to the bare quark mass. In the 
PNJL model the bare quark mass is a direct input parameter and the pion mass arises from the model so they are 
already related. The PNJL model yields a relation for the masses that is consistent with the lattice result (|9|. A 
comparison between these two is shown in Figure [T| where the PLSM line now corresponds to the lattice formula (|9| . 
The deviation between the models is small, but grows slowly towards higher niq being « 10 % at = 60 MeV and 
« 20 % at TUq = 250 MeV. Note that we obtain this consistency without tuning the coupling or the cutoff. Actually, 
we have checked numerically that as a function of both G and A, the values in Table |T] give the best agreement for 
the pion mass with the lattice curve. Altering G or A from these values will increase to^ for a given and thus 
adjusting these parameters would only move m^(TOg) further away from the corresponding lattice curve. 

Now we have a parametrization of both models in a way which allows us to study chiral symmetry breaking through 
one parameter, namely the quark mass rriq. In the PLSM model we have exchanged this parameter with the mass of 
the pion through a formula consistent with lattice result since in the case of PLSM model this is most convenient as 
the pion appears as an explicit degree of freedom. We now turn to investigation of explicit chiral symmetry breaking 
and its effect on the thermodynamics of these models. 




FIG. 1: Pion mass as a function of the bare quark mass in the PLSM and PNJL models. 



III. EFFECT ON THERMODYNAMICS 



Using the parametrization described in the previous section we can now study the relation between the pion mass 
and the critical temperatures of the models. There are a priori two transitions in the models: The chiral transition 
due to (approximate) restoration of chiral symmetry at finite temperature and the deconfinement transition encoded 
into the parameters of the potential for the polyakov field £. To both of these transitions one can assign their own 
critical temperature. Since over most of the parameter space the transition is a crossover, the definition of a critical 
temperature is vague. We have taken as definition of the critical temperature of the transition to be the temperature 
in which the temperature derivative of the corresponding field has its peak value. This definition however has some 
problems since a rapid increase in the field will yield a peak in the derivative regardless of the absolute value of the 
field at this point. To illustrated possible ambiguities, consider Figure |2] where the polyakov field £ and its temperature 
derivative are shown. At large ruq a double peak structure is visible in the derivative with a broad peak and a sharp 
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peak (a similar effect arises also when ruq is less than 5 MeV, but not shown in Figure^. When one looks at the 
values of the field 0, the sharp peaks correspond to the discontinuous jumps in the field at large while the broader 
peaks signify the location of the largest continuous rise in the value of the field. It is not clear which of these rapid 
changes should actually be used to determine the critical temperature. Some insight to the problem can be found by 
looking at the Figure [3] where the constituent mass and its derivative are shown. One immediately notices that at 
large niq the discontinuities in the constituent mass in Figure [3] correspond to the jumps in the Polyakov field (f> in 
Figure [2j This suggests that the sharper peaks in the Polyakov field derivative are merely an effect produced by the 
rapidly changing constituent mass and its interaction with the Polyakov field. Also since the overall behaviour of the 
Polyakov field at different values of rUg is not significantly altered by the appearence of these sharp peaks it seems 
that the broader peak is a better indicator for the critical temperature. This can be further justified by considering 
that deconfinement is characterised by £ ^ 1 and the broader peak is roughly at ^ 1/2 which is the point where the 
system is becoming dominantly deconfined in contrast to being confined when ^ ~ 0. From here on we will call the 
temperature associated with the broader peak as the (primary) Tc and temperature associated with the sharper peak 
the induced Tc- We observed a similar double peak structure also in [6J and it seems to be a feature of these effective 
models, but may not be true in two-flavor QCD. 
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FIG. 2: The Polyakov field (j> and its derivative in the PLSM model for various values of niq. The derivatives have been 
normalized to unity for clarity. 
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FIG. 3: Constituent mass and its derivative in the PLSM model for various values of niq. The derivatives have been normalized 
to unity for clarity. 

Note that the determined from the derivative of the chiral field is dependent on the bare quark mass, or 
equivalently on mj^. For PLSM this behavior can be affected by increasing the value of coupling g: larger the value 
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of g, slower the increase of Tc as m-^ increases. However, since g oa it turns out that our simple parametrization 
yields larger value for g than what one would get if one would fit Mn more carefully from lattice |16j . Projecting 
this on our analysis implies that while we observe PLSM and PNJL transition temperatures to behave qualitatively 
similarly up to at least to^ — 800 MeV, it is plausible that quantitatively the critical temperature in PLSM might 
grow more rapidly with than in PNJL model. Larger values of g also lead to sharper transitions and may be 
responsible for the discontinuous transitions encountered in the PLSM case. 

Turning to the PNJL case, then, we first observe that the double peak problem does not arise here, although 
observed also in the PNJL model at finite /x in A related difference with respect to the PLSM case is observed 
comparing Figures |3] and [5] The constituent mass behaves very differently in the two models. In the PNJL model the 
chiral transition softens as we go towards higher mg which can be seen as a broadening of the peaks in the derivative. 
In the PLSM case, on the other hand, the transition was observed to become sharper and finally discontinuous as 
ruq is ramped up. It can be attributed to this difference that the double peak structure in the Polyakov field arises 
in the PLSM model but not in the PNJL model. Despite the difference in the sharpness of the chiral transition the 
transition temperatures determined from the peaks of the derivatives are relatively close to each other. 




FIG. 4: The polyakov field (j> and its derivative in the PNJL model for various values of m,. The derivatives have been 
normalized to unity for clarity. 




FIG. 5: Constituent mass and its derivative in the PNJL model for various values of m,. The derivatives have been normalized 
to unity for clarity. 

Let us then sketch the phase diagrams of these theories. The phase boundaries of both models are shown in Figure 
[6] for a range of pion masses. From previous studies we know that when m^r = 138 MeV both, chiral and polyakov, 
transitions occur at a temperature of 211 MeV in the PLSM model and at 229 MeV in the PNJL model. Away from 
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the physical point these two models behave similarly as a function of pion mass to about m^r = 800 MeV from which 
upwards the PLSM chiral Tc has a rising curve while the PNJL model curve starts to level. Also the Polyakov field 
behaviour is similar in both models, although it seems to be more strongly coupled to the chiral sector in the PLSM 
case when m^r < 600 MeV. At larger m^r the Polyakov field 'decouples' from the chiral sector in both cases and its 
Tc seems to level off at near the pure gauge value. Note that in Figure |6] for the field (j) in the PLSM, the critical 
temperature is the primary Tc discussed above. The induced Tc follows the chiral transition curve and is not shown 
separately. In addition to the PLSM and PNJL models Figure [6] shows the critical temperatures for the bare NJL 
and LSM models as well as for the pure gauge potential. Comparing the NJL and pure gauge critical temperatures 
with the PNJL model, one can see that the critical temperatures are shifted closer together at small ttt-tt. One can 
also see that the interaction induces a dependence on the Polyakov field critical temperature which is independent 
of TTiTr in the pure gauge case. Also naturally the dependence of the chiral field is altered. At large 711^^ both of 
these effects die out as the different transitions separate. The chiral part of PNJL approaches the NJL curve and 
the Polyakov field draws near to the pure gauge curve. The same things can be said from the PLSM case where the 
chiral and Polyakov sector critical temperatures are shifted towards each other at small mjr, but at larger separate 
and approach the bare LSM and pure gauge curves respectively. Only difference being again the stronger interaction 
between the two sectors in the range 200 — 600 MeV of the pion mass. 
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FIG. 6: The critical temperatures of the chiral and the deconfinement transitions shown against the pion mass in both models. 



IV. CONCLUSIONS 



We have considered the PLSM and PNJL models at non-physical pion masses in order to determine how these 
models react to explicit chiral symmetry breaking. The motivations for this study were in general terms to determine 
how well these models encode the QCD dynamics of chiral symmetry breaking and, in more specific terms, to test if 
the intertwining of deconfinement and chiral symmetry restoration in these models survives away from the physical 
value of the pion mass. For the PLSM model this analysis was made possible by a lattice-based parametrization which 
allowed us to set the model parameters in a consistent way away from the physical pion mass. In the PNJL model 
the pion mass arises from the model and can be changed by adjusting the bare quark mass which is a direct input 
parameter of the model. As we have shown, both models lead to similar behavior for m^(m^) and are in this extent 
consistent with existing lattice data. 

In section |TTT] we studied the thermodynamics, or more precisely the transition temperatures and phase diagrams, of 
the two models as a function of pion mass. We found that the behaviour of the polyakov sector is almost identical in 
the two models, as expected. The only real difference being the notable interaction between the chiral and polyakov 
sectors at low pion masses in the PLSM model. However the overall behaviour of the polyakov field was not largely 
affected by this interaction and the primary Polyakov Tc agreed with the one acquired from the PNJL model. In 
contrast to the Polyakov sector, the chiral transition temperatures showed a different behaviour in the two models 
at pion masses larger than 800 MeV, but this should rather be taken as indication of the fact that such large pion 
masses are beyond the validity of simple chiral models. 

For the pion masses above 400 MeV our analysis for the PLSM model confirms the results found in [5J and 
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qualitatively the results obtained for the PNJL model are similar. Since such masses make the applicability of these 
chiral models questionable, one should concentrate more on the range below ttItt = 400 MeV. In this range we 
observe that in qualitative agreement both PLSM and PNJL models imply coincident chiral symmetry breaking and 
deconfinement as well as modestly increasing critical temperature as a function of the pion mass. Within this range 
the models agree quantitatively within ^ 10% although it should be noted that the critical temperature grows more 
rapidly with increasing pion mass in PLSM model than in PNJL model. Since the quark mass can be thought of a 
more basic parameter of these models than the pion mass and since the pion mass corresponding to given quark mass 
differ in these two models, it is informative to also plot Tj, as a function of the quark mass. This is shown in Figure 
[7] and clearly confirms the result that near the chiral limit the crossover transition Tc can be determined either using 
the chiral or Polyakov loop degrees of freedom. As one moves further away from the chiral limit, above ruq ~ 100 
MeV corresponding to 711^^ ~ 400 MeV (see Figure [ij, the description using chiral fields becomes more uncertain. The 
results obtained from the Polyakov loop for Tc are similar in the two models which is not surprising since the mean 
field potential for the Polyakov loop is the same in both models and only weakly affected by the value of the quark 
mass through the interactions with chiral fields. 




FIG. 7: The critical temperatures of the chiral and the deconfinement transitions shown against the quark mass in both models. 

Existing lattice data for two-flavor QCD thermodynamics at pion masses above 400 MeV [17^ together with our 
results in this mass range suggests that the dependence of the critical temperature on the pion mass is better described 
by the Polyakov loop dynamics. In the model framework considered here this is natural consequence of the underlying 
coupled dynamics. Lattice data closer towards the physical value of the pion mass would be highly desirable in order 
to better constrain these effective models. 
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